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Abstract 

The physical properties induced by a quenched surface magnetic field in the Ising 
model are investigated by means of boundary quantum field theory in replica space. 
Exact boundary scattering amplitudes are proposed and used to study the averaged 
quenched correlation functions. 
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1 Introduction 



There has recently been increasing interest in the exact estimations of boundary effects in 
statistical models and quantum field theories, following seminal works by Cardy based on 
conformal invariance as well as a pioneering paper by Ghoshal and Zamolodchikov ||, 
based on the integrability of the boundary interactionF]. New results have been obtained in 
many subjects, as for instance, the subjects of quantum impurity problems and dissipative 
quantum mechanics || [7], [8| or the exact calculation of correlation functions and Casimir 
energy in the presence of boundary conditions |TJ, [| |T(], [TT|. Aim of this paper is to enlarge 
the range of applicability of boundary field theory to systems which present boundary 
effects induced by a quenched surface disorder and to propose an exact non-perturbative 
approach based on boundary scattering theory]]. The model considered here is the two- 
dimensional Ising model in the euclidean half-plane T = {(x,y) G R; x > 0} with a 
random magnetic field coupled only to the boundary spins. This model was previously 
investigated by renormalization group techniques and conformal field theory methods in 



14] . Its dynamics may be described by the euclidean action 

r+oo 

A = A + dyh(y)a(y) , (1.1) 



where Aq is the bulk action relative to the Ising model in its high-temperature phase 
in a half-space geometry, with free boundary conditions adopted at x — 0. In eq. fll.lj ) 
<j(y) is the local boundary magnetic operator and h(y) is a quenched random field with 
Gaussian distribution of mean h and variance A: 



-2 



h(y) = h ; h( yi )h(y 2 ) = h + A 5( Vl - y 2 ) . (1.2) 

In the bulk, the massive excitation of the model may be regarded as an interacting 
bosonic particle A(6) created by the magnetization operator]]. Its elastic S- matrix is 



simply given by S = — 1 ||18|| . The case A = corresponds to the usual Ising model with 
boundary magnetic field, which was originally discussed by Ghoshal and Zamolodchikov 
][|: these authors proved that the model is integrable and therefore all its properties (as 
for instance, its partition functions and its correlators || [K|]) can be recovered by using 
both its exact 2-body 5*-matrix in the bulk and the elastic boundary scattering matrix 

1 See also H for boundary scattering theory and refs. [Q, ^| for a general introduction to boundary 
effects. 

2 We refer the reader to the references [l7|] for all details and notations relative to exact boundary 

and bulk S'-matrices. 

3 As usual, 9 is the rapidity variable which parameterises the dispersion relations (E,p) — 
(m cosh 9, m sinh 9) of the particle. 
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given by 

7 . (6 m\ i sinh 6 — k , 

11(9, h) = ztanh [-- — ) . , 1.3 

K J V 2 4/i smh + k v y 
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where k = 1 — tt . 

2m 

In terms of microscopic processes, the scattering event is nothing else but the result of 
an infinite number of interactions which take place on the boundary and which enter the 
self-energy of particle A, each of these interactions involving the insertion of the magnetic 
field h (see for instance Fig. l.a for a graph of order h 6 ). The integrability of the model 
means that all scattering processes at the boundary proceed without particle production 
or absorption, i.e. an observer will always see a particle arriving at the boundary and 
bouncing back. 

The graph technique can be also used to account for the interaction due to the random 
fields. The only new feature is that we have to average over the random variables h(y), an 
operation which is performed^ by using eqs. (O). This can be graphically represented by 
pairing by means of dotted line the small circles which represent the magnetic fields (see 
for instance Fig. l.b), each of which gives a factor A. These extra terms will certainly 
increase the intricacy of the self-energy of the particle propagator but they cannot spoil 
the integrable nature of the boundary magnetic interaction, i.e. the one-to-one nature 
of the scattering process of the particle which hits the boundary. The reason is that the 
graphs which result from the average are similar to those entering the pure system which 
are integrable. Hence, we conclude that the dynamics associated with the Ising model 
with quenched boundary magnetic field may also be described by an integrable model. 

This paper is organised as follows: in Section 2 the model is formulated in the replica 
space and the exact reflection S'-matrix is computed. In Section 3 we compute some 
quenched averaged correlation functions using the exact scattering theory. In Section 
4 we apply the Thermodynamical Bethe Ansatz to compute boundary entropies. Our 
conclusions are in Section 5. 



2 Boundary scattering in replica space 

The quenched averages of the disordered systems can be obtained by using the so-called 
replica trick which is based on the identity logZ = lim n ^o(Z n ~^)/ n i where \ogZ is the 
free energy of the model under investigation. In this way the original random problem is 
mapped onto a quantum field theory problem, described by an effective action involving 
n degrees of freedom in the limit n — > 0. Apart from subtleties in taking the n — > limit, 

4 It may be convenient to separate the contributions which originate from the average value h by 
redefining h(y) — > h(y) — h. 
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the advantage of this transformation is that it makes available to us powerful tools of 
quantum field theory and therefore opens up a way of obtaining non-perturbative results 
on the original random system. For the problem of interest in this paper, it is easy to 



see that by using Q1.2Q , the effective action is given by 

r+oo 



/-pOO 
dy£*.(v) (2.1) 
i -°° a 

/+oo n 
a<b 



The action ( |2.1|) describes n copies of the Ising model hinged together at the boundary 
(Fig. 2). In addition to the coupling to the (mean) boundary magnetic field, the new 
interaction in the action ( |2.1|) describes a process where a particle with replica index a 
arrives at the boundary and then changes to a particle with another label b. This ele- 
mentary process may be repeated along the boundary an arbitrary number of times until 
the particle is finally re-emitted. Since no production events are involved, the boundary 
action ( |2.1|) defines to a boundary integrable theory. In the following we will assume the 
validity of the above argument irrespective of the actual value of n, in particular also for 
the values n < 2. Consequently, there are two possible channels the elastic reflections 
7tab{0) can go through: (1) the first channel is the reflection process described by the 
amplitude P(9) where the particle A a (8) bounces back by keeping its original replica 
index a; (2) the second channel is described by the amplitude Q(9) where the particle 
A a {9) changes its replica index a — > h as result of the interaction (Fig. 3). Therefore we 
write 

TZ ab {9) = P(9)5 ab + Q{9){\ - 5 ab ) . (2.2) 

In terms of the geometrical representation of Fig. 2, the dynamics consists in the free 
propagation of the particle on any one of the n sheets, according to its replica label: 
when the particle hits the boundary it can either remain on the same plane (with a 
probability | P{9) | 2 ) or redirect its motion to another of the n — 1 planes (with a 
probability | Q{9) | 2 ). 

The two amplitudes P{9) and Q{9) satisfy the unitarity conditions 



P{9) P{-9) + (n - 1) Q(9) Q{-9) = 1 , 

P{9) Q{~9) + Q{9) P{-9) + (n - 2) Q{9) Q{-9) = 



(2.3) 



as well as the cross-unit arity equations 




-P(f + 9) 
Q(f + 9) 



(2.4) 
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In the cross-unitarity equations we have chosen S a b = 1 — 25 a b as the iS"— matrix relative 
to the collision processes of the replica in the bulk. The system of equations ( |2.3| ) and 
( [2.4|) admits two different classes of solutions depending whether n ^ 2 or n = 2. 

For n / 2, the above crossing equations and the second of the unitarity equations 
may be identically satisfied by the following ansatzf] 



P(d) = (n ~ 2) (- - I) f(6) 



Q(0) = f(0) , 



(2.5) 



with f(6) being a crossing symmetric function 

f(0) = f{m - 6) . 



(2.6) 



Plugging this ansatz into the first unitarity equation we obtain the additional equation 

1 1 



/(*)/(-*) 



for the function f{6). The simplest "minimal" solution of ( p.6|) and ( |2.7| ) is given byQ 



(2.7) 



in 



2) ( e_ n 

' \iw 2(n- 



F(6) 



where 



F(9) 



r ( i 



-2) 



2-ki 



v ( « 

1 U(n- 



■2) 1 2ni 



■p I n I 1 I 6 \ p / n 6* 
V4(n-2) 2 27r V \4(n-2) 27T? 

The function F(8) admits the following integral representation 



{2.1 



(2.9) 



F(0) = exp 



00 (it e 2 ("- 2 ) 6t 
1 I n~ s i n tt~ 

t COSh 7 27T 



(2.10) 



Note that the above minimal solution does not contain any free parameters. Apart from 
possible CDD factors, this seems to be a general property of the system of equations 
( |2.3| ) and ( |2.4| ) for generic n. For n = 2, however, the situation is different. This is due 
to a simple reason: in fact, by unfolding two semi- infinite planes the system becomes in 
this case an ordinary Ising model (in the infinite plane) but with a line of defect. Hence 

5 A detailed discussion on the solution of eqs. ( |2.3| ) and ( |2.4| ) for n^2 can be found in Appendix A. 
6 The most general solution is obtained by multiplying (2.8) by an arbitrary CDD factor, i.e. an 



arbitrary meromorphic function < t > (6') which satisfies both the equations <t > (6>)$(— 8) = 1 and $(0) 
$(wr-0). 
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for n = 2 the system of equations ( [2.3|) and ( |2.4|) admits extra solutions, more precisely 
those which belong to the one-parameter family of ^-matrices of the Ising model with a 
line of defect which were determined in [[HJ . In the present notation they may be written 

as 

sin x cosh 9 



P{9) 
Q{0) 



sinh 9 — i sin \ 

(2.11) 

cos x sinh 9 



sinh 9 — i sin x 
where 

SinX=_ 7XA£- (2 ' 12) 
1 + 4 

These additional solutions at n = 2 may shed light on the identification of which physical 
situation the minimal solution (|2.8| ) corresponds to. Let us compare, in fact, the minimal 
solutions P{9) and Q(9) in the limit n —>■ 2 with the reflection amplitudes of the defect 
line: in this limit P(0) and Q(9) go to and to —1 respectively and therefore, from ( [2.11 ) 



and (|2.12j ) we see that they correspond to an S-matrix of the defect line with an infinite 
value of the coupling constant A. This comparison suggests that the minimal solution 
( p.8|) , without any adjustable parameter, may correspond to n Ising models infinitely 
coupled at the boundary, i.e. to the strongly disordered case A — > oo. 

Note that for n — 1 (and h — 0) we cannot of course have any interaction between 
the replica and correspondingly the diagonal amplitude P(9) in this case reduces to 

P(0) = -icoOi(j-^ , (2.13) 

which coincides with the reflection matrix of the pure Ising model with free boundary 
conditions. This result indicates that the simplest way to introduce the dependence on 
the (mean) magnetic field h is to multiply the minimal solution (|2.8|) by the CDD factor 

H(9) = tanh 2 — . (2.14) 

For n->0we have a well-defined limit and the exact boundary S"-matrix is given in this 
case by 

f) _ hl rf±--Mrf-^- N i 

p(Q\ _ 2 V 2 2m J \2iTi) TT(n\ 

V 2 ~ 2-Ki J V 2m, 



(2.15) 
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In the opposite limit n — > oo , we find that Q(9) goes to zero whereas P{9) becomes 

P W = -nanh^ T --j— — , (2.16) 

i.e. a pure phase. 

It is interesting to study the structure of the boundary bound states by varying n 
(in the following we will only consider the case h = since the discussion of the pole 
structure induced by the CDD term H{6) may be found in 0): from eqs. (|2.8|) and ( j2.9|) 
we find that the positions of the poles are at: 



«l = (a + i)„ + T (— 



a" o; ■ m ( n 
9, = —21711 I , 

1 2 \n-2j 



(2.17) 



where k, I = 0, 1, 2, .... For n > 2 there are no poles in the physical strip < Im 9 < tt/2. 

4fc+l < n < 8fc+4 
2/fc+l — IL — 4/c+3 



Vice- versa, for n < 2, as long as n belongs to the interval < n < |It| (respectively 



< n < H+l)) we nnc l ^ na ^ i n the physical strip there is the unique pole 9' k (respectively 
6"). These poles are associated to boundary bound states and their existence appears 
compatible with the analysis of the renormalization group relative to the variable A 
[ PH , @- In fact, for n > 2 A is a (marginal) relevant variable, hence of feeble effect 
at the short distance scales near the boundary; for n < 2, A becomes on the contrary 
a (marginal) irrelevant variable and therefore it may have quite a strong effect at the 
short distance scales near the boundary: this may result in an attractive force which 



can be strong enough to produce bound states. Note that the above poles ([2.17 ) are 
generically simple poles but for specific values of n a pair of them may coincide, giving 
rise to a double order pole. The residue at the simple poles is related to the (square) of 
the boundary-particle couplings (Fig. 4) 

R »(Q) „ l egato (2 lg) 

So, for instance, looking at the pole at 9 = —i~ (^z^J with < n < 1, we have for the 
boundary-particle coupling in the diagonal channel P(9) 

Y ( 2n-3 \ 

A. = 2 ^W^T • < 2 ' 19 > 

Analogous expressions are found for the residues at the other poles. 



6 



3 Quenched Correlation Functions 



Quenched averaged correlation functions of the model may be computed by using the 
replica formalism and quantum field theory methods. Here the calculation for the first 
few is addressed, namely the one and two-point functions of the energy and magnetization 
operators. To this aim the following general formulas are needed 



(<p(x)) = hm((p a (x)) 

n— »0 



(<p(x)<p(y)) = lim ((p a (x)(p a (y)) ; (3.1) 



(<p(x)) (<p(y)) = lim(<p a (x)(p b (y)) \ a 



In our case these formulas have to be further specified to take into account the boundary 
effects. Let us see how this can be done. 

In the presence of boundary interactions, a convenient way to perform the calculation 
of the correlation functions is to employ a method based both on the boundary state 
wave function and the form factors relative to the local operators in the bulkn: assuming 
that the boundary is placed at the "time" t — and described by its wave function | B >, 
the correlation functions may be in fact expressed as 



(0 


T t [O x {x u t x ) 


• • • Op {Xp, tpj\ 


B) 


(0 


B) 



(0 1 (x l} t l )...Op(x p} t p )} = 1 11 1V "/nim ' ( 3 - 2 ) 



Since in this geometric setting the Hilbert space of the theory is the same as that in the 
bulk, the local operators O a can be completely characterized by their bulk form factors 

(Ai^i),---, I 0.(0,0) | A ip+1 (9 p+1 ),...,A iq (9 q )) , 

irrespectively of the presence of the boundary. Hence, once the boundary state wave 
function is known, the calculation of the correlation functions ( |3.2| ) is in principle just a 
matter of introducing the completeness relationships of the intermediate states between 
the various operators. In our model the boundary state wave function is explicitly given 
by 



B) = exp 



1 r ~ X d6K ab {6)A a {-6)A b {6) 



0) , (3.3) 



.2 

where K ah {6) = P(6)5 ab + Q(0)(1 - 6 ab ), with P(0) = P (f - 9) and Q(0) = Q (f - 6). 
With the above information, let us proceed to the calculation of some correlation func- 
tions. 



7 For the model under investigation in this paper, the form factors are all known, being the form 
factors of the thermal Ising model. Their expression may be found for instance in ]l9j. 
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The simplest one is the averaged one-point function of the energy operator: in the 
language of the replica this is given by 

e (t) = <i(^) = lim (0 | e a (p) \ B) , (3.4) 

n — >U 

where p = (x, t) and t is the distance from the boundary. Once averaged, the one-point 
function (e(t)) does not depend on x, as a consequence of the translation invariance along 
the boundary. In the bulk, the operator e a (x,t) couples only to the two-particle state 
with the same index a and its exact form factor is given by |19| 

(0 | e a (x,t) | A a {9 1 )A a {9 2 )) = 2ixm% sinh e ±zll (3.5) 
x exp \—mt (cosh Q\ + cosh 9 2 ) + imx (sinh 81 + sinh 8 2 )} . 

Hence there is only one term of the boundary state wave function which contributes in 
this case (see Fig. 5) and we have 

(e(t))= -im / d6smh6P(9)e- 2mtcoshe . (3.6) 
Jo 

In the above formula, P{9) is understood to be evaluated^ at n — 0. Various profiles of 
this correlation function for different values of the mean magnetic field h may be found 
in Fig. 6. They present the typical cross-over associated to a relevant boundary operator: 
for h = the correlator present a short-distance singularity (e(t)) ~ — 1/ (2t) whereas for 
h = 00 there is a swap of the overall sign (e(t)) ~ l/(2t). For finite values of h (however 
large), the corresponding curves follow the behavior of the curve at h — at very short 
distances while they tend to follow the behavior of the curve relative to h — 00 at large 
distances. 

Let us consider now 

G(pu P2 ) = (e(pi) e(p 2 )) = lim (0 | e a ( Pl )e a (p 2 ) | B) . (3.7) 

n— »U 

In terms of the replica, there are four possible graphs entering the above correlation 
functions (see Fig. 7) and correspondingly we have 

G{p u p 2 ) = I x + h + h + h ■ (3-8) 
The first of them (Fig. 7.a) is just the product of the vacuum expectation values, 

h = eo(*i)e (t 2 ) • (3-9) 

8 Unless explicitly stated, from now on both P{9) and Q(9) are meant to be the functions evaluated 
at n = 0. 
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Let us adopt in the following the notation r = ti + t 2 , t = £1 — t<i, x = xi — £2, 
r = V^ 2 + r2 ■ The second term (see Fig. 7.b) is the bulk energy-energy correlation 
function, which, in terms of the Bessel function K may be written as 



d 



d 



—K (mr) + —K (mr) — m (K (mr))' 



dx 



dr 



(3.10) 



The remaining two terms relative to the graphs of Figs. 7c. and 7.d involve the amplitude 
P(9) and their expressions may be cast in the form 



( d \ ( d 

2m I — K (mr) \ F{x, r) — 2mK (mr) I — 



F(x,r) 



(3.11) 



(mF(x, r)Y 



' d_ 

dx 



F(x,T) 



' d_ 
dr 



F(x,t) 



d(5 P{(3) exp[— mr cosh (3 + imx sinh (5\ 



where we have introduced the auxiliary function 

Altogether, the two-point function ( |3T7| ) can be finally expressed as 

1 2 



(3.12) 



eo(h)e (t 2 ) + 



d_ 

dx 





2 











K Q (mr) + mF(x, r) 







+ 



mK (mr) + -—Fix, r 
ox 



d_ 

2 



K {mr) 



(3.13) 



In the case of the energy operator, it is also interesting to calculate the following quenched 
averaged correlation function 



Ga(puP2) = (e(ii))(e(t 2 )) = hm (e a {x)e b (x)) \ a 



(3.14) 



There are only two graphs which contribute to this correlation function. The first one is 
the disconnected term given by the product of the vacuum expectation values, eq. ( |3.9| ). 
The second one is the graph drawn in Fig. 8 which involves the off-diagonal amplitude 
Q{9). The correlator Ga{pi,P2) can be expressed in this case as 

\ 2 



Ga{P\iP-i) 







Z(x,t) 



dx 

+ eo(£i)e (t 2 ) , 

where we have defined the auxiliary function 

1 



d_ 
df' 



Z(x, t) 



(mZ(x, t)Y 



+ 



Z(x,t) = o / d(3Q{(3) exp[— mt cosh/? + imx sinh (3} . 



(3.15) 



(3.16) 
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In the limit £ — > the composite operator E a b(p) = lim^ e a(p + £) e b(p) does not present 
ultraviolet singularities: in fact, the two fields e a and e& do not have short-distance 
divergences in the bulk and they are linked to each other only through the boundary 
interaction. Therefore the above formula (|3.15|) can be specialized to the case p\ = p 2 



and studied correspondingly as one-point function of the field E ab {t): its profile versus the 
distance from the boundary is plotted in Fig. 9 for different values of the mean magnetic 
field. It is interesting to note the rapid crossover which occurs in this correlation function 
for large values of h at very short distance scales. 

Let us now turn our attention to some correlation functions of the magnetization 
operators. The simplest is the one-point function of the disorder operator p(p) 



p (t) = (p(p)} = lim(0 | p a (p) | B) . (3.17) 

In the high-temperature phase, the disorder operator p a {p) has a non-zero vacuum ex- 
pectation value and couples to all states with an even number of particles with the same 
replica index: its explicit form factors are given by[] ]19| 



(0\p a (0,0)\ A^)... A a (9 2n )) = (-^Utanh^-^ . (3.18) 

i<j 1 

Since the boundary state consists of a condensate of Cooper pairs, i.e. couples of particles 
with equal and opposite momentum, we have to specialize the above formula to the case 
(0 | /i a (0) | A a {— 9 1 )A a (9 1 ) . . . A a (—6 n )A a {6 n )) . This matrix element can be conveniently 
written as 



(0 | ^(0,0) I A a (—9i)A a (9i), . . . A a (-8 n )A a (6 n )) = % n (jltanh^J x det W(6 h 9 3 ) , 

(3.19) 

where W(6i, 8j) is an n x n matrix whose elements are given by 



(2 wcosh Qi cosh^,- \ 

The one-point function of the disorder operator is made of an infinite number of terms, as 
the one drawn in Fig. 10. Its final expression can be conveniently written as a Fredholm 
determinant 

<M*)> = £-r / de x ...de nZ n m ? tanh ^ p (^) e " 2mtcosh9fc det = 

n=0 711 J -°° \k=0 J 

= det(l + zW) , (3.21) 



9 The magnetization field denned by the Form Factors normalized as in eq. ( 3.18 ) differs from the 
corresponding conformal field only by a normalization constant, n(p) = T (J-conf(p), where T = 
2-i/i2 e i/8^-3/2 m -i/8 and A _ !_282427 (Glasher constant). 
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where the kernel of the integral operator is given by 

3 cosh6*i + cosh^ 

and 

E(9,mt) = e - mtcoshd (iP(0) tanh9) 1/2 , z = — . (3.23) 

As consequence of the translation invariance along the boundary, po(t) depends only on 
its distance t from the boundary. In terms of the eigenvalues of the integral operator and 
their multiplicity, p{t) can be also expressed as 



M*) = U(l + z^r • (3-24) 



As far as mt is finite, the kernel is square integrable and therefore its properties are 



those of a bounded symmetric integral operator For large values of mt } (n(t)) falls 
off exponentially to its bulk vacuum expectation value. However, when mt —>■ 0, the 
integral operator becomes unbounded: in this case, the multiplicity of the eigenvalues 
grows logarithmically as a(t) ~ ^ In and the one-point function presents a power 
law behavior (fi{t)) ~ A/(2t)^. To determine (, observe that as long as n ^ 2 (which we 
always assume to be the case in the following) and h ^ oo, we have lim^oo P(9) = —i. 
Let us denote this limit as P_. For n^2 and h = oo, we have lim^oo P(Q) = i, instead. 
This limit value will be denoted by P + . To study the short distance behavior mt — > 0, 
we can reasonably substitute P{9) with its asymptotic limits P±. The eigenvalues of the 
integral operator then become dense in the interval (0, oo) according to the distribution 

Hp) = — r — 

cosh ivp 

and for the exponent (± relative to the two cases we find 



C±= -nJo dph V±^J =-- + ^arccos( T l 



3 
8 
1 

8 



(3.25) 



The profile of this one-point function is drawn in Fig. 11 for several values of h: the 
observed cross-over behavior of the curves by varying h is perfectly analogous to the one 
of the energy operator. 

Let us conclude this section with the discussion relative to the (averaged) two-point 
function of the magnetization operator a(p). In the bulk, this operator has non-vanishing 
matrix elements on an odd number of particles and their expression is given by the same 
formula as (|3.18|) Consider first 



Q'fapi) = MpiMpz)) = lim (0 I a a ( Pl )a a (p 2 ) | B) . (3.26) 
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There are two kinds of graphs entering the above correlator (see Fig. 12). The first does 
not involve the boundary and therefore the sum of such graphs gives rise to the known 



correlation function in the bulk which is expressible in terms of Painleve' function [22|. 
The second set of graphs is made of an arbitrary number of pairs of particles A a emitted 
by the boundary and absorbed by the two operators. The two operators, in turn, are also 
linked together by an arbitrary (odd) number of intermediate particle states. The first 
non-trivial contribution to this correlation function is the one shown in Fig. 13, given by 

J( Pl ,p 2 ) = 2g 2 F(x,T) , (3.27) 

where g denotes the (constant) one-particle form factor of the magnetization operator 
(0 | cr(0) a | A b (3)) = g5 a b and the function F(x,t) is defined in ( |3.12| ). A convenient form 
of the series resulting from the graphs of Fig. 12 is presently unknown. 

Let us consider now the other type of two-point averaged correlation function 



Ga(PuP2) = (o-(Pi)) <o"(P2)> = lim (da(piK(p 2 )> \a^b • (3.28) 

n— >0 

The graphs contributing to this correlation function are those shown in Fig. 14. This time 
the boundary can emit an arbitrary number of pairs of particles of type A a A a , A b A b or 
A a Ab. the former two pairs are individually absorbed by each operator whereas the latter 
link the operators to each other. The simplest (non-disconnected) term of the resulting 
series is given in this case by 

J%\ Pl ,P2) = 2g 2 Z(x,t) , (3.29) 

where the function Z(x,t) is defined in eq. ( |3.16| ). The above discussion can be easily 
generalized to the correlator Ga(pu Pz) = (p(Pi)) (p(P2)) as we ^- m this case the lowest 
contribution of the series originates from a graph like the one of Fig. 8 and its explicit 
expression is given by 

1®( \ f + °° df]l f + °° d ^ 2 f VV A-/^ Afa x Afa w 

Ja (Pi,P2) = — tanh ( — ) Q{Pi) Qifo) x. 

JO Zll J-oo Z7T 2 

e _ m(tl+t2 ) (cosh ^ 1+cosh ^) co ^ m ( Xi _ ^^(sinh^ + s inh/3 2 )] . (3.30) 
It can be rewritten as 

4 2) ( Pl ,p 2 ) = 4 J2 [(2/+ 1) (z 2l+2 {x,f)Z 2l (x,f) -Z 2 2l+1 (x,f)) + (3.31) 



1=0 



(3.32) 



+ (21 + 2) (Z 2l+3 (x, t)Z 21+1 (x, t) - Z 2l+2 (x, t] 



where 



1 r +OD dB ( 3\ l - 

Zi(x,t) = - — tanh— Q(3) exp[— mt cosh 3 + imx sinh/3] . (3.33) 

2 J-oo 2n \ 2 J 
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Since the rate of convergence of the functions Zi(x,t) becomes faster and faster by in- 



creasing Z, the series (|3.31|) can be efficiently estimate by means of the first few terms. 



As was the case for the energy operator, we can define the operators 

E ab (p) = \ima a (p + £)a b (p) ; T ab (p) = \imp a (p + £)p b (p), (3.34) 

without facing any short-distance singularities. It is therefore possible to study the 
correlator G/L^iPii Pz) m the limit p\ — > P2, the so-called Edwards-Anderson (EA) order 
parameters (a(p)) 2 and (p(p)) 2 . The profile of the Edwards-Anderson order parameter 



(/i(p)) 2 for different values of h by using the lowest order approximation (|3.30| ) is drawn 



in Fig. 15. Presently however it is quite difficult to make any comparison with the lowest 
order expression in A (h — 0) for such parameter 



(Mt)> 2 ~ A-^ , (3.35) 

which was obtained by Cardy ||15|| . This difficulty is due to several reasons: (a) the 
different methods which were employed in the two cases (Cardy obtained in fact ( [3.35D 
by using boundary conformal perturbation theory; on the other hand, our approach is 
based on the scattering theory relative to the off-critical excitations); (b) the lack of an 
exact re-summation of the series originated from the graphs of Fig. 14 which, if known, 
would permit in principle a comparison between the two methods above and finally (c) 
that the result ( |3.35|) applies for small values of A whereas our scattering theory seems 
to apply to the opposite limit A — > oo instead. 



4 Boundary Thermo dynamical Bethe Ansatz 

Once the exact ^-matrix of a model is known, finite-size effects and associated thermo- 
dynamical quantities can be calculated by means of the Thermodynamical Bethe Ansatz 
(TBA): the case of systems with relativistic invariance and a cylinder geometry with 
periodic boundary conditions has been put forward in ||12|| . What is of interest for us 
here is the generalization of the TBA to systems with boundary: such generalization can 



be found discussed in 11 



Consider then an Ising model defined on a cylinder of width L and length R with 
boundary conditions at its extremities described by the reflection scattering theory of 
Sect. 2 (see Fig. 16. a). Choosing as direction of the time the horizontal axis between the 
two boundaries, the partition function of the model can be expressed as 

Z = (B | exp(-RH) | B) , (4.1) 
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where H is the Hamiltonian of the bulk system with periodic boundary conditions. In 
the limit mR 3> 1, the partition function reduces to 



Z~ gi g Iie - H ^ , (4.2) 

where E is the ground-state energy of the system whereas gj and gjj are the boundary 
degeneracy at each end of the cylinder. The /^-independent term 



hi<7/,// (4.3) 



can be then interpreted as boundary entropies of the system [[n| . They can be computed 



by using the boundary Thermodynamical Bethe Ansatz fly]. The application of this 
method requires first of all to look at cylinder geometry of Fig. 16. a in a different way, 
namely the time evolution should be regarded as it takes place along the vertical axes, 
i.e. along the circumference of length L. Hence, the dynamics consists in a set of TV 
particles which move and scatter each other along the segment of length R until they 
reach the boundaries: then, they are reflected off with amplitudes lZ a b{6). The finite 
geometry of the system induces a quantization condition on the momenta of the particles 
which can be obtained as follows. Let n a be the number of particles of replica index a, 
with J2 a n a — N. Consider now one of N particles, of species b and send it on a "round 
trip" along the cylinder (Fig. 16. b): after coming back to the same location, the wave 
function of the N particles has picked up a new phase shift. In fact, each time that the 
particle scatters off a like particle it picks up a —1 from their bulk S-matrix (all other 
scattering processes in the bulk with particles of different species do not contribute since 
the S'-matrix in these cases is simply +1); in addition to these processes, there is an 
additional phase Rbk(0)Rkb(0) (sum on the internal index k) due to the scattering off the 
two boundaries. Since the total wave function of the iV particle state is assumed to be 
periodic under this "round trip" along the cylinder, we get the quantization condition 

e 2iRm b smh(6)f^\n b +n k -lf^(Q\\2 _ ^ (4 4) 

where A a (9), a = 1,2, ... ,n, are eigenvalues of the reflection matrix TZbk(9). They can 
be found explicitly since any n x n matrix of the form 

My = (P- Q)6 tJ + Q 

has only two types of eigenvalues: Ai = P+{n— 1)Q with multiplicity one and A 2 = P—Q, 
which is (n — 1) degenerate. 

With this information, the analysis that follows is quite standard (see for instance 



and therefore we only sketch the main steps here. Eq. (4^4) puts a constraints 



among the physical rapidities and therefore for each particle in the interval 9 and 9 + d6 
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we can define the density of holes, Pi(9), and of actually occupied states, Pi{6). Hence, 
the density of available states for particles of species a is p a {9) + p^{9) which is given by 

2n(p a (9) + p h M) = m a cosh(^) + ^a{9) , (4.5) 
where we have introduced 

$„(0) = -2i^ In A a (9) - 2ttS(9) . (4.6) 

The delta function term is introduced in order to remove the unwanted solution 9 = 0. 
The partition function is obtained by expressing the Helmholtz's free energy as a func- 
tional of the density of states (holes and particles) and minimize it. The final expression 
of the boundary entropy is given by the massless limit of 

S b = lim V / -rMO) Ml + e~ eM ) , (4-7) 

1=1 

where <S& is the boundary entropy at one end and the pseudo-energies £i(8) are given by 
the simple expression 

€i{9) = miRcosh(9) . (4.8) 

since the bulk S-matrix is ±1. 

The formula (14.71 ) becomes then 

n r°° I 

S h = - dx-— 2 ln(l + e~ 2 ™ x ) , (4.9) 

7T JO 1 + X 1 
2 

where mRcosh(9) = 2nax and a = h i?/(47r). Note that each species gives exactly the 
same contribution and so we have an overall factor of n in our expression. The above 
integral can be computed exactly 



*«- (;)')■ w 

For n = 1 (which corresponds to the pure Ising model) we obtain the correct boundary 
entropy difference ASb = lny^ (for a = oo and a = 0) |T5[ . The entropy difference of 
the disordered model (n = 0) is computed by taking the derivative of the ground state 
degeneracy with respect to n and let n — > 0. This gives ASb = In a/2, exactly the same 
result as in the pure Ising model. 

5 Conclusions 

In this paper we have proposed an exact scattering theory in the replica space of n species 
of particles A a (in the limit n — > 0) to describe the dynamics of the two-dimensional Ising 
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model coupled to a boundary random magnetic field. By using methods borrowed from 
boundary quantum field theories, we have then computed averaged correlation functions 
of the order parameters as well as finite-size effects of the system defined on a strip. It 
would be obviously interesting to test the feasibility of this approach for analyzing other 
models with random interactions localized on the boundary. 
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Related work, on the 2-dimensional Dirac theory in semi-infinite space with ran- 
dom boundary interactions, has recently been done independently by S. Guruswamy and 



A.W.W. Ludwig p 
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Appendix A 



In this appendix we briefly discuss the most general solution of the boundary unitarity 
and cross-unitarity equations ( |2.3|) and (|2.4j) for n^2. Let's initially define a function 
D{9) = P{6) — Q{9). It's easy to see that the unitarity equations impose 

D(9)D(-9) = 1 . (A.l 

Therefore, substituting P{9) = D(9) + Q(9) into the first unitarity equation (|2.3| ), we 
find 

D{9)Q{-9) + D(-9)Q(9) + nQ{9)Q{-9) = . (A.2 
Dividing by Q(9)Q(—9), we have 

D{9) , D{- 



Q{9) Q(- 

Define T{9) by 



+ n = 0. (A.3 



From ( |A.3|) we see that T{9) is an odd function. Cross-unitarity requires 

T{m-9) +T(9) =n-2 . (A.5 
With the position -D(#) = d(9){2T(9) — n), we see that d{9) has to satisfy 

= (2 m - „m(-e ) -„ ) ■ (A - 6 

From this point on the analysis is standard and one can easily find two solutions for d{9 
(up to CDD factors) 

r/rfll = 1 4(n-2)_M2 4(n-2) j / A ? 

w (2T(g) - n) rq - n+2T ^)ya _ "- 2T W) ' 1 

v v 7 ' v 4(n-2) / V2 4(n-2) / 
, m - 1 + 4(n-2) Ml 2 + 4 (n-2) ^ , 

U (2rw + n)r(i + ^)r(i + ^)- ( ' 

These are general expressions that solve unitarity and cross-unitarity conditions provided 
T{9) is an odd function satisfying the equation ( |A.5|) . The simplest solution for T{9 
which gives the correct analytic behavior is 

T(9) = (n-2)—. (A.9 

in 

With this choice of T(9), the two solutions equations ( |A.7| ) and ( |A.8| ) are easily seen to 
be connected by a CDD factor. They differ in the positions of the poles: the first one has 
poles in the physical strip for n > 2, the second one, instead, for n < 2. This last feature 
is one of the reasons to use the latter solution for the physical scattering amplitudes in 
the replica space. 
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Figure Caption 



Figure 1 . Graphs of the perturbative series: (a) pure system (b) system with disorder. 

Figure 2 . Replica space. 

Figure 3 . Reflection amplitudes in replica space. 

Figure 4 . Boundary bound state and boundary-particle couplings. 

Figure 5 . Graph of the averaged one-point function of the energy operator. 

Figure 6 . Profiles of the one-point function of the energy operator versus its distance 
from the boundary for different values of k: k — 1 (long dashes), k — — 1 (short 
dashes) and k — — 10 (full line). 

Figure 7 . All possible graphs of the averaged two-point function of the energy operator. 

Figure 8 . Diagram relative to (e(pi) e(p 2 ))- 

Figure 9 . Profiles of the one-point function (e(t)} 2 for different values of k: k — 1 
(long dashes), k — —1 (short dashes) and k — — 10 (full line). 

Figure 10 . One of the diagrams entering the one-point function of the disorder oper- 
ator. 

Figure 11 . Profiles of the one-point function (fi(t)) at the lowest order approximation 
for different values of k: k = 1 (long dashes), k = — 1 (short dashes) and k = — 10 
(full line). 

Figure 12 . Diagrams entering the averaged two-point function of the magnetization 
operator. 

Figure 13 . Simplest graph contributing to the averaged two-point function of the 
magnetization operator. 

Figure 14 . Diagrams contributing to the two-point function (fi(pi)) (fj,(p 2 )). 

Figure 15 . Profiles of the EA parameter (p(t)) 2 at the lowest order approximation for 
different values of k: k = 1 (long dashes), K — —1 (short dashes) and k = — 10 (full 
line) . 

Figure 16 . Strip geometry (a) and "round trip" (b) of the TBA quantization equation. 
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